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Abstract. This is the third and last part of a series of 3 papers. Using the same method 
and the same coordinates as in parts 1 and 2, rotating dust solutions of Einstein's equations 
are investigated that possess 3-dimensional symmetry groups, under the assumption that 
each of the Killing vectors is linearly independent of velocity and rotation w°' at every 
point of the spacetime region under consideration. The Killing fields are found and the 
Killing equations are solved for the components of the metric tensor in every case that 
arises. No progress was made with the Einstein equations in any of the cases, and no 
previously known solutions were identified. A brief overview of literature on solutions with 
rotating sources is given. 

I. Summary of the method. 

This paper is the third and last part of a series of 3 papers, for parts 1 and 2 see Refs. 
1 and 2 (Paper 1 and Paper 2). For convenience of the readers, this section is repeated 
after Paper 2. 

This is a concise summary of results that will be used in this paper. For proofs, 
motivations and references see Paper 1^. 

Every timelike vector field of unit length that has zero acceleration and nonzero 
rotation defines the functions r(x),?7(x) and ^(x) such that: 

Ua = r,a+'n^,a- (1-1) 

These functions are defined up to the transformations: 

T = T'-s{e,v'), ^ = F{e,v'), v = G{e,v'), (i-2) 

where the functions F and G obey: 

F,^/ G,n' —F,r,' G,^i = 1, (1.3) 
(this guarantees that the Jacobian of the transformation is 1), and S is determined by: 

S,^, = GF,^, -tI, S„, = GF„, . (1.4) 
If u°' is the velocity field of a fluid whose number of particles is conserved: 

(v^nn"),, = 0, (1.5) 

(where g is the determinant of the metric tensor and n is the particle number density), 
then one more function ({x) exists such that: 
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V^nu"^e"^^'^,fsVnC:S: (1-6) 
and it is determined up to the transformations: 

C = C' + T(e',V)- (1-7) 

Note that n is not defined uniquely by (1.5). For example, if = 6"o and n obeys 
(1.5), then n' = nf{x,y,z) (where / is an arbitrary function) will also obey (1.5). This 

nonuniqueness allows for a greater freedom in the choice of C than (1.7), and the freedom 
will be used in some cases. 

The following relations hold: 

0{t. //.(.() 



The last of (1.8) guarantees that {T,^,r),C} can be chosen as coordinates, they will be 
called Plebahski coordinates. Then, with {r, ^, 77, (} — x^, x^, x^} — {t, x, y, z}: 

900 = 1, fi'oi = y, 902 = 903 = 0, g = det{gai3) = -n~^, 

= nSl Uafs = -oo^a = il/2)5\5^p, (1.9) 

where is the rotation vector field, and ujap is the rotation tensor corresponding to the 
velocity field Ua' 

i^aP = ^{U'a,(3 - U(3,a " UaU/S + Uf^Ua), = -{'^ / V-g)^'^'^^^ UpUJ^s. (1.10) 

If ujai3 7^ and it" = (what is assumed throughout), then necessarily the pressure p — 
const and may be interpreted as the cosmological constant {k := SttG/c^). 

If any Killing vector field exists on a manifold (on which all the assumptions specified 
so far are fulfilled), then, in the coordinates of (1.9), it must be of the form: 

A;" = (C + - + (f>,y - 0„ 5^2 + XS^3, (1.11) 

where C is an arbitrary constant and (j){x,y) and X{x,y) are arbitrary functions of two 
coordinates. Whenever 0^q, 7^ 0, a transformation of the class (1.2) - (1.4) can be found 
that leads to: 

fc° = (5°i. (1.12) 

The metric then becomes independent of x, and the coordinates preserving (1.12) are 
determined up to the transformations: 
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t' = t- lyH,ydy + A, x'^x + H{y), y' = y, z' = z + T{y), (1.13) 



where A is an arbitrary constant and H, T are arbitrary functions. 

The condition „ 7^ that allows one to fulfil (1.12) means that the Killing vector A;" is 
linearly independent of the vectors and w" at every point of the spacetime region under 
consideration. In Paper 1, solutions of the Killing equations and of the Einstein equations 
were considered under the assumption that there exist three Killing vector fields on the 
manifold, two of which have = const in (1.11), while the third one has 7^ and can 
be transformed to the form (1.12). In Paper 2, it was assumed that only one Killing field 
has = const, while two have 0^q, 7^ 0. In the present paper, all three Killing fields will 
be assumed to have (p^a 7^ 0. One of them (/i;(i)) can be transformed to the simple form 
(1.12), while the remaining two will have the general form (1.11). 

In this Paper 3, no progress was made with the Einstein equations in any of the cases. 
Also, no related results were found in the literature except in case 1.1.2.2, see at the end 
of sec. V. 



II. The Lie algebra of the symmetry group. 

According to the assumptions made in the preceding section, there exist the following 
three Killing vector fields: 



= (C2 + - y<P,y + (P,y - ^2 + H^, y)^^ . 

kfs) = {Cs + ,p- y^j.y + 0„ 5^ - 6^ + A3(a;, y)Sl (2.1) 

where C2 and C3 are arbitrary constants, and 0, ip, A2 and A3 are unknown functions of 
{x,y), to be determined from the commutation relations. The coordinates of (2.1) are 
determined up to (1.13). 

The fields A;(i), k(^2) ^-nd A;(3) will form a Lie algebra if constants a, . . . ,j exist such that: 

[^(1), ^(2)] = ak^i) + hki2) + c/c(3), 
[^(1), ^(3)] = dk(x) + ek(2) + /A;(3), 

[k{2) , k^3)] = 9k(i) + hk{2) + j A;(3) , (2.2) 
Eqs. (2.2) are equivalent to the following set: 

yV,2/), (2.3a) 
(2.3&) 
(2.3c) 
(2.3rf) 

yip,y), (2.3e) 



0,x -y(i>,xy = &(C2 + - #,2/ ) + C(C3 + -0 - 

(p,xy^ a + b(p,y+cip,y , 

b(f),x+C1p,x , 
^2,x = b\2 + CA3, 

ip,x -yi>,xy = e(C2 + - y(l),y ) + /(C3 + i>- 
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i^,xy^ d + e(t),y+fil),y, (2.3/) 

i^,xx = e0,^ +fi^,x , (2.39-) 

A3,. - eA2 + /A3, (2.3/1) 

= /l(C2 + 0-#,J +j(C3 + ^- 1/^,2/), (2.3i) 

0,2/ V',X2/ - V',2/2/ - ^^,2/ 0,^2/ + "0,^ 0,2/2/ = ^ + ^0,2/ +^"0,2/ , (2-3j) 

-0,2/ + 0,a; V',X2/ + "0,2/ 0,xx " 0,a;2/ = -^0,x , (2.3A;) 

0,y A3,a; - 0,a; A3,j, - ^',2/ ^2,x + ^2,5/ = + jAs. (2.3Z) 

The equations (2.3a-c) are integrated with the result: 

0,^= ay + 60 + cV' + 6C2 + cC3; (2.4a) 
the equations (2.3e-g) are integrated with the result: 

= dy + ect> + fi, + eC2 + /C3; (2.46) 
and the equations (2.3i-k) are integrated with the result: 

(p.y -(f),x i>,y = gy + h(f) + j'tp + hC2 + jCs. (2.4c) 



The equations arc now sorted as follows. Eqs. (2.3d) and (2.3h) form a set that determines 
A2 and A3, cqs. (2.4a-b) form a set that determines and ip, the remaining two equations 
((2.4c) and (2.31)) are consistency conditions to be imposed on the solutions of the two 

former sets. 

Note that the set (2.4a-b) and the set {(2.3d), (2.3h)} are of the same general form: 
each of them is an (ordinary differential) linear vector equation of first order: 

U,x^AU + W, (2.5) 
where, for (2.4a-b), the constant matrix A and the vectors U and W are: 

while the set {(2.3d), (2.3h)} is homogeneous, so 1^ = in (2.5), the matrix A is the same 

as in (2.6) and C/ = |^ j . 

With the constants 6, c, e and / being all arbitrary, several cases will have to be 
considered separately. Just as in Paper 2, the cases that arise will be organized into a 
binary tree an numbered in a positional system that will enable one to quickly identify the 
complementary part of each alternative (see also the diagram). 

The first alternative appears in solving the characteristic equation for the matrix A. 
Its eigenvalues are: 
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where: 



(2.7) 



A := (6 - /)' + 4ce, e^^ 1,62^-1. (2.8) 

We first consider: 

Case 1: A 7^ (i.e. A has two distinct eigenvalues). 

The second ahernative appears immediately in finding the eigenvectors of A: the cases 
c 7^ and c — have to be considered separately. 

Case 1.1: c ^ 0. 

The solution of the set {(2.3d), (2.3h)} is then: 

A2-2cL2(z/)e'^l^ + 2cL3(z/)e"^^ 

A3 = (/ - 6 + VA)L2{y)e^^^ + if - b - VA)L,{y)e^-^ , (2.9) 

where L2{y) and L^{y) are arbitrary functions. The cases A > and A < could be 
considered together for a large part of the reasoning. When A < 0, cti and 0^2 are complex 
and a2 — oT. Then, L2 and L3 have to be complex, too, with L3 = L2. However, the two 
cases lead to different sets of Bianchi types, and so it will be convenient to split them here. 

Case 1.1.1: A > (i.e. both eigenvalues aie real). 

Then, in solving the set (2.4a-b), the cases detA 7^ and det74 = have to be 
considered separately. In the end, however, the case det A = Q turns out to be empty, i.e. 
in all subcases that arise in it there exists a linear combination of the Killing vectors A;(2) 
and /c(3) with constant coefficients that is spanned on u and w. This means that all these 
subcases are in the domain of Paper 2 and need not be considered here. Therefore we 
will do away with the case det A — Qhj only indicating the method of verification of the 
statement above. 

When det^ = 0, the following is true (from (2.7) - (2.8)): 

e = 6//c, A = (fe + /)^ = & + a, q;2 = 0. (2.10) 

Since we are still in Case 1 in which a2 7^ 0.1 by assumption, we can take it for granted 
that q; 7^ here. The solutions of (2.3d), (2.3h) and (2.4a-b) are: 

A2 = 2c[L2{y) + L3(y)e--], A3 = 2/L2(y) - 2bLMe"\ 

(f) = F(y)e"^ + (a/ - cd)xy/a + P{y) - {ab + cd)y/a'' - C2, 

^ = (//c)F(y)e"^ - {b/c){af - cd)xy/a - {b/c)P{y) 

-{f/c){ab + cd)y/a^-C^, (2.11) 

where F{ii) and P{y) are other arbitrary functions. The further procedure goes exactly as 
for the case det A 7^ presented below, and leads to the results specified above. 
Prom now on, in Case 1.1.1 we assume that: 
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det A ^ bf - ce ^ 0. 

Then the solution of (2.4a-b) is: 



(2.12) 



= 2cF(|/)e"i^ + 2cP(|/)e"^^ - ^^^y - C2, 

oj - ce 

= (/ - 6 + x/A)F(y)e"^- + (/ - 6 - \/A)P(?/)e"^^ - ^^^^y - C3. (2.13) 

OJ - ce 

The solutions (2.9) and (2.13) must now obey the consistency conditions (2.4c) and (2.3Z). 
Eq. (2.4c) is a polynomial in 6°^^ and e"^^ whose coefficients are functions of y, the 
polynomial contains e("i+"2)a; gcni;^ gQ!2a; terms independent of x. The cases a2 7^ — cti 
and a2 = —ai require separate consideration. We first consider: 

Case 1.1.1.1: a2 7^ —ai. 

Comparison of coefficients of e("i+"2)a; both sides of (2.4c) leads to: 

-a2F,yP + aiFP,y^0. (2.14) 

If F = 0, then (2.14) is fulfilled identically, and this case has to be considered separately. 
The result is similar as in the case det A = 0. The consideration, parallel to the one that 
follows below, reveals that with F = either one of the Killing fields becomes coUinear 
with rotation (and this situation is in the domain of Paper 2) or the symmetry group 
becomes two dimensional (which case is not considered here at all). Consequently, we 
proceed assuming F 7^ 0. The solution of (2.14) is then: 

P = /?F"2/"i, (2.15) 

where f3 is an arbitrary constant; ai 7^ because det A ^ 0. 
The coefficients of Fe"'^^ on both sides of (2.4c) imply: 



h = Mibf - ce)][-l(a/ -cd){f-b + VA)-ae + bd]-^{f-b + VA). (2.16) 
zc Ac 

In considering the coefficients of e"^^ in (2.4c) we have to set aside the case P = for 
separate consideration because the terms with e"^^ all vanish identically when P = 0. 
However, the case P = is in fact empty in the same sense as the case P = 0: either 
the Killing field /c(3) becomes coUinear with rotation, and this situation is in the domain 
of Paper 2, or the symmetry group becomes two-dimensional. Hence, we shall follow the 
case P 7^ only. Then the coefficients of Pe"^^ on both sides of (2.4c) imply: 

J = {bf - ce)-^[-{af - cd)f + (-ae + bd)c], (2.17) 
and the terms independent of x imply: 

g^{bf - ce)-^[-{af - cd)d + (-ae + bd)a]. (2.18) 
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Eq. (2.3/) is a polynomial in exponential functions of x that involves e^"^^^, e^"^^^, e("i+"2)a;^ 
e"i^ and e"^^. In consequence of (2.9) - (2.18) most of the resulting equations are fulfilled 
identically, the only one that brings in new information is: 

-2a2F,y L3 + 2a^FL^,y + 2(/3a2F"^/"0(^2i^,y /F - L2,y) = 0. (2.19) 
This is integrated with the result: 

L3 = (1/q;i)F"2/"i(7 + Pa2L2/F), (2.20) 

where 7 is a new arbitrary constant. 

Finally, the following Killing fields resulted: 

1^(1) — "1 ) 



k%) = [2c{F - yF,y )e"i^ + 2c/?F°2/"i(l - ^yF,y / F)e''^^]5^ 

+ {2cF,y e"i" + 2c/3— e"^^ - ^^r^)5^ 
ai of — ce 

-(2caiFe"i^ + 2ca2/5F°2/°^e"2^)5^ + pcLse"^^ + 2— ^"2/^^1(7 + /3a2i>2/i^)e"^^]53 , 



^"3) = [(/-& + VA)(i^ - )e"^^ + {f-b- v^)/3F-^/"Hl - /F)e"^-](5o" 

oil 

+ [(/ - 6 + \/A)F,y e'^'^ + {f-b- VA)/?— F"^/"i-iF,^ e"^^ - ^^^^]<^f 

-[(/ - 6 + v^)Q;iFe°i^ + (/-&- V^)/3Q;2F"2/''^e°^^]5^ 

+ [(/-&+ \/A)L2e"i" + {f-b- \/A)— F"2/"i(7 + /3a2^2/i^)e"^1(53". (2.21) 

These formulae are simplified by changing the basis in the algebra of the Killing vectors. 
Taking k[^^ = kf^^ + fjE^^^i) and k[^^ = kf^-^ + =^^k'^^^ instead of k'^^^ and kf^^ we obtain 
the same result as if: 



a = d = 0. (2.22) 

With (2.22), we take k'^-^ = {2^/A)-^[kf^) - ^(/ - 6 - \/A)A;^2)] instead of kf^-^ and k[q^ = 
(2c/3)~^(/c"2) — 2c^'"3-)) instead of /c^g)- Further simphfication results from the transformation 
(1.13) with: 

H^a^HnF, T,y= L2/{aiF). (2.23) 
The Killing field k^^^ does not change, while the other two become (all primes dropped): 

ki2) = e"''^{So - 0^26^ + [7/(/^«i)]^3 }, ^"3) = ^"'"(^0 - «i^2 )• (2.24) 



Using (2.22) in (2.16) - (2.18) we obtain: 



g = h = j = 0. (2.25) 

Looking at the commutation relations one can see that in effect we have also achieved 
c — e = 0, ai = b, a2 = f, even though the initial basis (2.21) was calculated under 
the assumption c 7^ 0. This allows us to predict that in the case c = 0, set aside for 
separate consideration, we will obtain (2.24) again; see case 1.2 in sec. VI. In view of the 
assumptions made earlier, we have: 

b + f^O^b-f, b^O^f, (2.26) 

and consequently the Bianchi type is VI^, with the free parameter being gb = {b + f)/{b — 
/)■ 

The Killing fields are simphfied even further after the following transformation that 
leads out of the Plebanski class: 

t' = {b- fr'e-f^bt + y), x' = X, y' = {b- fT'be-'^ft + y), 

z'^-^[(3{b-f)]-\t + y/b) + z, (2.27) 
that results in (primes dropped): 

fcfi) = -m + ^?- byS^, = 5^, = 5^, 

u"^{b- /)-^[6e-^^(5^ + bfe-^'^S^ - {-f/P)S^], = nS^. (2.28) 

In the new coordinates, the Killing equations imply for the metric: 

^00 = e'f^[l - U/bf + i^g,^/{b(5f - 2f^h23/m + /'M> 
^01 = ef^[^g^s/m - fhul 902 = e^'^f^^i^h^^/m - (6 - f)/b^ - //i22], 
^03 = e^^[7W(^'/9)- A23], ^12 = e''^/ii2, g22^e^'^h22, to = e^^/i23, (2.29) 
where g^n, /ii2, giz, /i22, ^23 and ^^33 are arbitrary functions of z. 

III. Ceise 1.1.1.2: a2 = — ai- 

All equations up to (2.13) still apply, but (2.4c) has to be reconsidered. With — — cci 
the following equations hold: 

/ = A = 4(62 + ce), ai = v^/2 = -q;2. (3.1) 

In considering eqs. (2.4c) and (2.3/) it may be assumed that F 7^ 7^ P because the 
opposite cases lead, just as before, out of the domain of this paper. With F 7^ 7^ P, 
the coefficients of e"i^, of e""^^ and the terms independent of x (which now include the 
coefficients of e("i+"2)a;j jj^ (2.4c) imply, respectively: 

h = A-^/2[-(a6/c + rf)(-26 + a/A) - 2(-ae + bd)] - j{-2b + VA)/(2c), 
j=a, cAFP={l/2){g + A-^[-h{ab + cd) + a{-ae + bd)]}y^ + B, (3.2) 
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where B is an arbitrary constant. 

In eq. (2.3/) only the terms independent of x provide a piece of new information, the 
other parts of the equation are fulfilled identically in consequence of (3.1) - (3.2). The new 
information is: 

Ls = ^/{aiF) - PL2/F (3.3) 

(the integration constant 7 was chosen so as to correspond to (2.20)). Like in sec. II, we 
first construct the Killing fields by substituting (3.1) - (3.3), (2.9) and (2.13) into (2.1), 
then simplify the result by changing the basis in the Lie algebra and carrying out the 
coordinate transformations (1.13) (the difference is that here the coefficient in the formula 
for k'^^ is A instead of (2c/3)~^). The result is: 

a = d = 0, A;fi) = S^, kf^^ = e^'^iS^ - a,S^), 
ki2) = e-'^'^ii-gy'' + 2B)5^ + 2gy5^ + a,{gy'' + 2B)5^ + 2{(rfA/a,)5^]. (3.4) 
The commutation relations are: 

[%),^(2)] = -Qii^(2), [k{2),k^3)]^2gaik(^i), [k^^), k(^i)] ^ -aik(^3), (3.5) 

and they correspond to the Bianchi type VIII when g ^ and type VIq when g — (note 
that tti 7^ by the assumption defining case 1.1.1). 

The Killing fields are further simplified by the coordinate transformation: 

t' ^t + y/ai-[2Bai/(c-fA)]z, x' ^ x, y' ^ y, z' 
that results (with primes dropped) in: 

^fi) = ^2) = e-"^^[2yy5r + ai(yy^ + 2B)b'^ + bf^, 

= 5^, = [nai/(2c7A)](-4S5o" + bf), 

goo = 1, 501 = y, 5-02 = -l/ai, 5-03 = 4S, (3.7) 

the other components of gij being just unknown functions. The Killing equations for A;"^^ 
now imply that the metric is independent of x, while those for A;^^ are solved by: 

gu = {aiyfg22 - 2aiyHi2 + 

g\2 = -oi\yg22 + -^"12, gi^ = -(yiyg23 + H13, (3.8) 

where Hu, H12, H13, g22, 923 and 5(33 are arbitrary functions of t and z. 

In solving the Killing equations for /c^^, three cases have to be considered separately: 

Case I: gB ^ 0. 

In fact, the cases gB > and gB < lead to different results, but the formulae for 
gB < can be easily reconstructed from those for gB > by taking real combinations of 



= [ai/{2c^A)]z (3.6) 
k%) = -a.e'^'^S^, 
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the complex solutions. Hence, we shall only present the formulae for gB > 0. The solution 
of the Killing equations is: 

/J, :— {2gBY^'^ai, U := /112 cos(4/x2;) — /ci2 sin(4/x2;), 

V := hi3 cos(2/iz) — ki3 sm{2fiz) 

Hii = -{2Bq\/ij)U + hii, H12 = /ii2 sin(4/i2;) + ki2Cos{Afiz) - Hi^/{ABai), 

Hi3 = hi3 sm{2iiz) + ki3 cos{2fj,z), §22 = {g/lAU - [g/{'^Baiij)]V + /i22, 

923 = \}Ji/{2Bal)]V - h33/{ABai), ^33 = h^, (3.9) 

where the hij{t) and kij{t) are arbitrary functions. 

When gB < 0, is imaginary. Then the trigonometric functions go over into the 
appropriate hyperbolic functions, and hu and hi3 have to be taken imaginary, too. 

Case II: B = 0. 

The solution of the Killing equations for /c^^ is then: 

Hu = hssiaizf + 2aihi3Z + hu, 
H12 = -ajghssz^ - Saighi^z'^ + (01/123 - 2ghu)z + hi2, i^is = ^ih^^z + his, 
922 = {oiig)^h33Z^ + iaig^hnz^ - 2g{aih23 - 2ghii)z^ - ighuz + h22, 

g23 = -aigh33Z^ - 2ghi3Z + /i23, 5'33 = ^33, (3.10) 
where the hij{t) are arbitrary functions. 

Case III: g = 0. 

The Killing equations for /c^^ imply here: 

Hu = [{4:Baifh22 + SBaih23 + h33]{aizf + {8Baihi2 + 2hi3)aiz + hu, 
H12 = {4:Baih22 + h23)oiiz + /ii2, H13 = (4Sq;i/i23 + h33)aiz + hi3, 

922 = h22, g23 = h23, g33 = ^.33, (3.11) 

where the hij{t) are arbitrary functions. 

With this, case 1.1.1 is exhausted and we go back to (2.9) to consider the other branch 
of the alternative. 

IV. Case 1.1.2: A < (i.e. the eigenvalues ai and a2 are complex and 
conjugate to each other). 

As stated before, it is more convenient to reparamctrizc (2.9) so that it contains only 
real quantities, and then repeat the procedure of case 1.1.1 in this new parametrization. 
We denote: 

y/^^D, F^G + iJ, P^G-iJ, L2^M + iN, Ls^M-iN, (4.1) 
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where G, J,M and N are new unknown functions of y, and then: 

a,,2^l{b + f) + Us,,2D, si = l, S2 = -l. (4.2) 
In this notation, eqs. (2.9) and (2.13) adapted to the case A < are: 

A2 = 4ce(''+^)^/2[Mcos(Da;/2) - Nsm{Dx/2)], 

A3 = ^(/ - b)X2 - 2De^^+f^''/^[Msm(Dx/2) + N cos(Dx/2)], 
2c 

(f) = Ace^^+f>/^[Gcos(Dx/2) - Jsm(Dx/2)] - - C2, 

oj - ce 

^ = 2(/ - 6)e(^+^)^/'[G'cos(L>x/2) - Jsin(L>x/2)] 

_2i;e(''+/)^/2[G'sin(L>a;/2) + Jcos(Da;/2)] - ~7 ^^"^ y - C3. (4.3) 

OJ - ce 

Just as in Case 1.1.1, these expressions must now satisfy the consistency conditions (2.4c) 
and (2.3/). However, in considering them, the cases b + f ^ and b + f — have to be 

taken separately. 

Case 1.1.2.1: 6 + / ^ 0. 

Both sides of (2.4c) are then polynomials in e^^'^^^^^'^, and some of their coefficients 
involve cos(Da;/2) and sin(Da;/2). The coefficient of e^''^-^^^ leads to the equation: 

AcD[-D{GG,y +JJ,y ) + {b + f) {GJ,y -G,y j) ] = 0. (4.4) 

We are working in the case c ^ ^ D, so only the expression in square brackets can 
vanish. Its form suggests the substitution: 

G = KcosL, J = K sin L, (4.5) 

where K and L are new functions of y. Eq. (4.4) becomes then —DKK,y +{b + f)K^L,y — 
0, and its solution is: 

K = Be^'^+fWD^ (46) 

where B is an arbitrary constant; B or else we are back in the domain of Paper 2. Eqs. 
(4.5) and (4.6) provide a parametric representation of G and J in terms of the function L, 
which is arbitrary at this stage. 

The coefficients of e^^^^^l"^ in (2.4c) involve sin and cos that always go in fixed pairs. 
The coefficients of {e(''+-^)^/2[G'cos(Da;/2) - Jsin(Da;/2)]} imply: 

h = [4c(6/ - ce)]-\af - cd){b^ -f + D^) + 2c(6 + /)(-ae + bd)] -j{f - 6)/(2c), (4.7) 
and the coefficients of {e^^+f^''/^[G sm{Dx/2) + Jcos{Dx/2)]} imply: 

j = -{bf - ce)-^[f{af - cd) - c{-ae + bd)]. (4.8) 
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Finally, the terms independent of x imply: 

g = -{hf - ce)-^[h{af - cd)+j{-ae + bd)]. (4.9) 
With (4.7) - (4.8), eq. (2.3/) is reduced to: 

4cl>e(^+-^)^ [-D{GM,y +G,y M + JN,y +J,y N) 

+ + f)iGN,y -G,y N - JM,y +J,y M)] = 0. (4.10) 

Only the expression in square brackets can vanish in the case now considered. After (4.5) 
and (4.6) are substituted into (4.10), the resulting equation integrates to: 

e-(''+^)-^/^{-M[L>cosL+(6 + /)sinL] + 7V[(6 + /)cosL-L>sinL]} = 7^(6 + /) = const. 

(4.11) 

Since h + f ^ Q ^ D hj assumption, this can be solved for N: 

N =[{h + f) cos L - D sin L]-^{^B{h + f)e-^^+fWD ^ L^{p^f^ L]}. (4.12) 

The resulting KiUing fields are: 

^{1) — "i > 

= 4ce(''+^)^/^[(G - yC,, ) cos{Dx/2) - (J - yj„ ) sm{Dx/2)]5^ 
+{4ce(^+/W2[G,^cos(Da;/2) - J,^ sin(Da:/2)] - ff^}^? 

OJ C€ 

-2ce^^+f^''/^{{b + f)[G cos{Dx/2) - Jsin{Dx/2)] - D[Gsin{Dx/2) + J cos{Dx/2)]}5^ 

+4ce(^+/W2[Mcos(Da;/2) - N sm{Dx/2)]S^, 

^^3) = 2e(^+^)^/2{(/ - b) [{G - yG,y ) cos(Dx/2) - (J - yj,y ) sin(Da:/2)] 
-D[{G - yG,y ) sin(L'a;/2) + ( J - yj,y ) cos{Dx/2)]}S^ 
^2e(''+/W2{(/ _ 6)[G',^cos(L>x/2) - J,ySm{Dx/2)] 

-D[G,ySm{Dx/2) + J„ cos(Da;/2)]} - "T^^'^ K 

OJ C6 

-e^^+M^{{f -b^- D^)[Gcos{Dx/2) - Jsm{Dx/2)] 
-2fD[Gsm{Dx/2) + J cos{Dx/2)]}S^ 
^g(6+/)x/2|2(j _ 6)[Mcos(L'x/2) - Nsm{Dx/2)] 

-2D[Msm{Dx/2) + N cos{Dx/2)]}6^. (4.13) 
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Again, the formulae simplify when the basis in the Lie algebra is changed. First, we take 
k'»^ = kf^^ + f^/e^i) and k'^-^ = /c^g) + =ifz^k'^i) instead of A;^2) and /c^g^ respectively. The 
result is equivalent to: 

a = d = 0, (4.14) 

and then (4.7) - (4.9) simplify to: 

g^h^j^O. (4.15) 
With (4.14) taken into account, we change the basis again by taking k'^-^ — {—2BD)^^{k^^-^ — 

^/c°2)) and A;(2) = {4:cB)~^k'^2) instead of k^^-^ and k^2), and carry out the transformation 
(L13) with: 

H = 2L/D, T,y = 2[B{b + f)]-^e-^^+f^^^^{McosL + NsinL) (4.16) 
The result is equivalent to L = = M, i.e.: 

G^B, J = 0, N^-fB. (4.17) 
The resulting KiUing fields are: 

t.a _ ra 
/C(l) — Oi , 

A;^3) = e(^+^)^/'[sin(L>x/2)5^ - ^VS^ + -fC0s{Dx/2)S^, (4.18) 

where: 

W := (6 + /)cos(Da;/2) -Dsin(Dx/2), 

V := Dcos(Dx/2) + (b + f) sm(Dx/2). (4.19) 
The commutation relations are: 

[h,k2] = ^{b + f)k2-^Dks, [k2,h]=0, [k3,k,] = -^Dk2-^{b + f)k3, (4.20) 

and they correspond to Bianchi type VII/i with the free parameter as — —{b + f)/D. 
Since [k2, k^] = 0, coordinates can be adapted to /c2 and ks simultaneously. The following 
transformation does it: 

t' = -D-^e-^^+f^''/^ [Wt + 2 cos{Dx/2)y] , x' = x, 
y' = e-(^+-'^)^/2|^-iy^ _ 2W-^[1 - D-^ cos{Dx/2)V]y}, 

z' = {^/D){b + f)t + 2{^/D)y + z, (4.21) 

but the new coordinates are no longer in the Plebahski class. After the transformation, 
with primes dropped: 
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kil) = li-ib + f )t + Dy]6^ + 6^ - ^[Dt + {b + f )y]6^, kf,^ = 6^, k^,^ = 6^, 

u"" ^ D-'e-^''+^^''/^-W5^ + VS^) + {-f/D){b + f)5^, w'' ^ nS^ (4.22) 

The formulae for goa (a = 0, 1, 2, 3) in terms of gij = 1, 2, 3) in the new coordinates 
are given in Appendix A. 

In the new coordinates, the metric is independent of t and y, while the Killing equations 
for /c"^^ have the following solution: 

511 = hn, 912 = e(''+^)^/'[iy/ii2 - (l/D)(b + f) cos{Dx/2)his], = h,^, 



g,, = e(^+^)-{[f (6 + fn,,/D' + 1] cos\Dx/2) 

-2{-f/D){b + f) COs{Dx/2)Wh23 + W^h22}, 



g,, = e(^+^)^/^[iy/i23 - {l/D){b + f) cos(L>x/2)/i33], 533 = ^^33, (4-23) 

where the hij{z) are arbitrary functions. The components goa in terms of those given above 
are given in Appendix A. 



V. Case 1.1.2.2: / = -b. 

We go back to eqs. (4.3) which simplify as follows: 

A2 = 4c[Mcos(Da;/2) - iV sin(Da;/2)], 
A3 = -(6/c)A2 - 2D[Msm{Dx/2) + N cos{Dx/2)], (5.1) 

= 4c[G'cos(L'a;/2) - Jsin(L'a;/2)] - Z y - C2, 

b^ + ce 

iP = -46[G' cos(Da;/2) - Jsin(Da;/2)] 

-2D[G'sin(L'x/2) + Jcos(L'x/2)] + ~f ^ ^"^ y - C3, (5.2) 

0^ + ce 

With these new forms of 0, ■?/',A2 and A3 we reconsider (2.4c) and (2.3Z). Eq. (2.4c) is 
now linear inhomogeneous in sin(Da;/2) and cos{Dx/2). The coefficients of cos{Dx/2), of 
sin(Dx/2) and the terms independent of x imply, respectively: 



ab + cd 2 



D'G - 2aDJ = 4chG - AbjG - 2jDJ, (5.3) 



b^ + ce 

D'J - 2aDG = -AchJ + AbjJ - 2jDG, (5.4) 



ab + cd 2 



6^ + ce 
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-4cD^{GG,y +JJ,y )^gy+{b'^ + ce)-^ [-h{ah + cd) + j{-ae + hd)]y. (5.5) 

From (5.3) and (5.4) it follows that 2(j - a)D{G'^ + J^) = 0. Since D hy assumption, 
and G = J = leads to the domain of Paper 2, this implies: 

j = a. (5.6) 
With (5.6), eqs. (5.3) and (5.4) reduce to: 

h - [4c(62 + ce)]-\ab + cd)D^ + ab/c. (5.7) 

The integral of (5.5) is: 

-2cD^(G^ + J^) ^{g + (6^ + cey^[-h(ab + cd) + a(-ae + bd)}y^/2 + B, (5.8) 
where B = const. 

Eq. (2.3/) leads now to two additional equations, one of which has the solution: 

GM + JN^-f/D^ const, (5.9) 
and what remains of (2.3Z) is then: 

4aD[M sm{Dx/2) + N cos{Dx/2)] = 0. (5.10) 

This has two solutions, a — and M — N — 0, that must be considered separately. 
However, the case M — N — turns out to be included as the subcase 7 = of the 
formulae below. Hence: 

a = j = 0, h=^{b^ + ce)-'^dD'^. (5.11) 
The resulting KiUing fields are: 

A;f2) = 4c[(G' - yG,y ) cos{Dx/2) - (J - y J„ ) sm{Dx/2)]S^ 

cd 

+{Ac[G,yCos{Dx/2) - J,ySm{Dx/2)] - -^^K 
+2cD[Gsm{Dx/2) + J cos(Dx/2)]5^ + Ac[M cos{Dx/2) - N sm(Dx/2)]5^, 

kfs) = {-46[(G - yG,y ) cos(Dx/2) - (J - yj,y ) sin(L>x/2)] 
-2D[{G - yG,y ) sm{Dx/2) + ( J - yJ,y ) cos{Dx/2)]}6^ 
+{-4:b[G,yCos{Dx/2) - J,ySm{Dx/2)] 
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-2D[G,yShi{Dx/2) + J,yCos{Dx/2)] + -g^^}^! 

+{2bD[-Giim{Dx/2) - J cos{Dx/2)] + D'^[G cos{Dx/2) - J sm{Dx/2)]}6^ 

+{-4b[Mcos{Dx/2) - Nsm{Dx/2)] - 2D[M sm{Dx/2) + N cos{Dx/2)]}6^. (5.12) 

By changing the basis to /c^") — ^(2) + (^^ + ce)~^cdA;"i) and k'^-^ — k^^-^ — {b^ + ce)~^bdk^i>^ 
the result equivalent to: 

d^h^O (5.13) 

is achieved, and then, from (5.8): 

+ J2 = -(gy^ + 2B)/(AcD^). (5.14) 
Eq. (5.14) suggests the parametrization: 

G^KcosL, J^KsinL, (5.15) 

then, from (5.14): 

= G^ + J^ = -{gy^ + 2B)/{AcD^), (5.16) 

and L remains arbitrary. With a = d = we change the basis again to k"'(2) ~ (4c)~^/C(2) 
and k"fs) = {~'^D)^^[k'^^ + {b/c)k'^^] and we carry out the transformation (1.13) with: 

H = 2L/D, T,y= -2{MsmL- NcosL)/{DK). (5.17) 
The result is equivalent to L = = which implies: 

J = 0, G = K, M = -f/{DK), (5.18) 
and the Killing fields become: 

'^{1) — "1 > 

A;") = {K - yK,y ) cos{Dx/2)6^ + K,y cos{Dx/2)6^ 
+^DKsm{Dx/2)5^ + [-f/{DK)] cos{Dx/2)5^, 

A;«^ = (X - yK,y ) sm{Dx/2)5^ + K,y sm{Dx/2)S^ 

~DKcos{Dx/2)S^ + [y/iDK)] sm{Dx/2)S^. (5.19) 
The commutation relations are: 

[fci,fc2] = -^^3, [^2,^3] = -b/(8cD)]fci, [fc3,fci] = -^m2 (5.20) 

16 



The Bianchi type depends on the constant g: for g^/c > it is type IX, for (7/c < it is 
type VIII and for 51 — it is type VIIq. The last case is contained in (4.20) - (4.23) as the 
subcase 6 + / = 0. 

The Killing equations for k"^^ imply that the metric tensor is independent of x. Knowing 
this, one can simplify the Killing equations for fc^^ and k'^^y Since /c^^ and /c^^ are linear 
in sin(Dx/2) and cos(Da;/2), while g^p are independent of x, each Killing equation implies 
two equations: the coefficients of cos{Dx/2) and of sin(Da;/2) have to vanish separately. 
The pair of equations implied by k^2) is identical to the pair implied by fc^-j, and it is: 

(K - yK,y )g^p^t + b/{DK)]g^p,, + ^D^KS'^g^/s 

+5j{-yK,yy g^o + K,yy gal - bK,y /{DK^)]gas} = 0, (5.21) 

Kga/3,y + Sl{-{K - yK,y )gop - K,y gip - [-^ /{DK)]g^p} 
+K,y 5lg2p + 5l{-{K - yK,y )g«o - K,y g^^ - [7/(L'X)]g«3} + K,y 5}ga2 = 0, (5.22) 
The solution of (5.22) is: 

gn = + 2Hss{'yK/Df J K-^R{y)dy + 2{'y/D)K'HisR{y) + K^Hi^, 
gi2 = {j/D)H23R{y) + H,^, gn = {j/D)H,sKR{y) + KH^s, 

g22 = H22/K^, g23 = H23/K, g33 = H33, (5.23) 

where the Hij{t, z) are arbitrary functions, and R{y) is: 

R{y) := j K-'dy. (5.24) 

With K given by (5.16), R{y) and / K^^R{y)dy can be easily calculated, but the result 
has to be given separately for gB ^ 0, ioi g = and for 5 = 0, so eq. (5.24) is the most 
compact notation (but see below). 

Note that 7 and B cannot vanish simultaneously; if 7 = = S, then K — yK,y = = 
K^yy, and then eq. (5.21) implies gi2 = g22 = 923 = 0. Together with gQ2 — (we are 
still in the Plebahski class) this means that det{gai3) = 0. Also, g and B cannot vanish 
simultaneously because with g = = B we are back in the domain of Paper 2. With 
^2 _|_ ^2 _^ 0, the following new variables can be introduced for solving (5.21): 

u = 2cD-ft + Bz, -Bt + 2cD-fz. (5.25) 

With 7 and B running through all possible values, the hypersurfaces u — const are timelike, 
null or spacelike. However, the solution of the Killing equations has the same dependence 
on u and v in every case. 

For solving (5.21), the cases gB ^ and gB = have to be separated. 
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Case I: gB ^ 0. 

(this means we are considering the Bianchi types IX and VIII here, but not VIIq) . In this 
case: 

R= -2cD^y/{BK), (5.26) 
and the solution of (5.21) (with (5.23) - (5.24) aheady taken into account) is: 

:= gB/{85^D^), 5"^ := {B/Df + {2c^f, U := hu sinh(2A^;) + ku cosh(2A^;), 

Hu = -cD^U/{25^X) + /ill, Hu = hi2 cosh(2At;) + sinh(2Ai;), 

iJi3 = -[cD^/{25^\)][h2^sinh{\v) + ^23 cosh (At;)], 

H22 = -2S'XU/{cD') - [gB/{2c'D')]hn - [Sc^V iBD^)]hss, 

H23 = /i23 cosh(Av) + /C23 sinh(Av), i^aa = /ias, (5-27) 

where the hij{u) and kij{v) are arbitrary functions. 

Eqs. (5.27) are adapted to the case gB > 0. When gB < 0, A^ < 0, i.e. A IS imaginary. 
Then ku and ^23 have to be taken imaginary (and U becomes imaginary in consequence 
of this), the hypcrbohc functions go over into the corresponding trigonometric functions 
in the well-known way, and the functions if^ remain real. 

Case II: gB = 0. 

Then K,yy = 0. The form of R{y) still depends on whether g ^ 01 g — and 5 7^ or 
B = 0. The formulae below apply in each case. The solution of (5.21) (again with (5.23) - 
(5.24) already taken into account) is here: 



icygyY , 4 c^^Dgy^ o /C^-D^, cjgy , . o cL>^ , , 

rr ngy u 2 cD^ . . TT {lay? ,2 lay ^ , , 

^13 = ^^I^hs3V - J^Kh23V + /lia, H22 = ^^^^^hsV - -^h2sV + h22, 

H23 = - ^^^3^ /^33^ + Kh23, = hs. (5.28) 

The subcase of (5.28) in which g = and the hypcrsurfaces u = const are spacelike should 
have a common subset witht the class considered by Demiahski and Grishchuk^. These 
authors considered Bianchi type VIIq models with nonzero rotation, with spacelike orbits 
of the symmetry group which are flat and with the source being a perfect fluid (the pressure 
is not constant in their class) . A member of the present collection should result when p — 
const. However, Rcf. 3 docs not contain sufficient information to identify it. 

With this. Case 1.1 is exhausted. We go back to (2.8) with A and consider: 



VI. Case 1.2 : c = 0. 

Just as it was announced in the paragraph after (2.25), this case brings no new infor- 
mation. Three situations occur here: 
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1. The group becomes two-dimensional (because two Killing vectors become coUinear); 
these cases are not considered here. 

2. A linear combination of the Killing fields with constant coefficients is spanned on 
u°' and w°'; these cases are in the domain of Paper 2. 

3. In the case when the group is 3-dimensional and none of the Killing fields is spanned 
on u"' and w", the formulae are equivalent to (2.25) - (2.29) and (3.2) - (3.11) (both sets 
reappear) . 

The proof consists simply in retracing the whole reasoning from (2.8) on with c = 0. 
As seen from (2.8), with c = necessarily A > 0, and so no analog of Case 1.1.2 arises 
here. The essential steps of the reasoning are described in Appendix B. 

Case 1 is exhausted at this point. 



VII. Case 2: A = (i.e. A has one double eigenvalue). 

The reasoning from (2.8) on has to be repeated with this new assumption. Then the 
double eigenvalue is a and: 

6 = /±2v^, a = /±v^=(6 + /)/2. (7.1) 

The double sign denotes two different cases, but they will be considered at one go. In 
finding the solutions of (2.5), the cases deiA ^ and det^ = have to be considered 
separately. We first consider: 

Case 2.1: detA ^ 0. 

This means: 

h + f^Q^a. (7.2) 

In the next step, the case c = has to be set aside for separate consideration, so we first 
assume: 

Case 2.1.1: c 0. 

The solutions for A2,3, and are here: 

A2 = {[1 + (6 - f)x/2\L^{y) + cxL,{y)}e^'^^^^'\ 
As = {-{h - ffxL-,{y)/{Ac) + [1 + (/ - 6)x/2]L3(t/)}e(''+^W^ (7.3) 



= {[1 + (6 - f)x/2]F{y) + cxP{y)}e^'+f>'' + A{cd - af)y/{b + ff - C^, 

^ = {-{b - ffxF{y)/{Ac) + [1 + (/ - h)x/2\P{y)}e^'+f>l^ 

-[a{h - ff/c + Ahd]y/{h + ff - Cs, (7.4) 

With such A2,3, and t/j, eq. (2.4c) involves polynomials of second degree in e^^'^f^^^'^ and 
also of second degree in x. The equations imphed by the coefficients of x^e^^+-^^^ and of 
2;g(b+/)a: are fulfilled identically, while the remaining ones are: 

-(Ac)-\b-f fFF,y+fF,yP- bFP,y- CPP,y=0, (7.5) 
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[{b - f)F/2 + cP][d-h+{a + j){h - /)/(2c)] = 0, (7.6) 

[{h - f)F/2 + cP]{[a{h - 3/) + Acd]/{h + f) - j} = 0, (7.7) 

g = {-Ahicd - af) + Mb - /)Vc + 4M]}/(6 + ff. (7.8) 
The vanishing of the first factor in (7.6) and (7.7) leads to the relation: 

hs) + ^^(2) c(6 + /) - + ^'^^ ' 

which means that the combination on the left is collinear with w", and so this case belongs 
to the domain of Paper 2. Hence, (7.6) and (7.7) imply: 

h^d+{a + j){b- /)/(2c), j = [a{b - 3/) + Acd]/{b + /). (7.9) 
The solution of (7.5) may be represented parametrically by: 

F ^2(AR + cR\nR)/(b + f), P ^ -(b - f)F/(2c) + R, (7.10) 

where A is an arbitrary constant and R{y) is an arbitrary function. It can be assumed 

that R ^ because otherwise we are back in the domain of Paper 2. 

Eq. (2.3/), with the functions given by (7.3) - (7.4) is a polynomial of the same form 
as (2.4c), and only the coefficients of e^'''^-^^^ provide a new equation: 

-(4c)-^(6 - ff{F,y L, + FL2,y) + fF,y L3 - - c{L,,yP + LsP,y ) 

-bL2P,y+FL2,yP = 0. (7.11) 

Using (7.10), the solution is found again in a parametric form: 

L2^2S{A+c]nR + c)/{b + f) + ^R, = -{b - f)L2/{2c) + S, (7.12) 

where 7 is another arbitrary constant and S is another arbitrary function. The resulting 
Killing fields arc now calculated from (2.1) using (7.10) and (7.12). Since this procedure 
has already been performed a few times in this paper, we shall not quote the intermediate 
results. It turns out that ^(3) contain terms which arc constant multiples of k"^^ , 

these are removed when a new basis, k'^^ and k'^^^ is appropriately defined. Then we change 
/c^g) again, to /^''^g^ = k'^-^ + [{b — f) / {2c)]k'^^ and carry out the transformation (1.13) with: 

H^2\nR/{b + f), T^2{b + f)-'l {S/R)dy. (7.13) 
The Killing fields that result are (with primes dropped): 

Ua _ 
/C(l) — Oi , 

k^2) = e(''+^)-/2{[cx + 2A/{b + f)]S^ -[c + A+{b + f)cx/2]S^ + ^S^} 

kf,^^e^'^f^^/'[5^-^{b + m]. (7.14) 
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It can be assumed now that .4 = bccaTisc this is equivalent to changing the basis to 
^(2) — ^{2) ~ 2^/c^)/(6 + /). The commutation relations are: 

[hi),k{2)] = ^{b + f)k(^2) + ck(3), = 0, [k{3),k(i)] = -^{b + f)k(3), (7.15) 

and they correspond to Bianchi type IV. 

Because of [A';(2), A;(3)] = 0, coordinates can be adapted to /c(2) and /c(3) simultaneously. 
The following transformation does it: 

t' = e-(''+-^)^/2 + 1 (5 + + a-^] ^ ^> ^ 
2 

y' = -c-'e-^'+f^^/'i^ib + f)t + y], z'^ 7(& + f)t/{2c) + ^y/c + z, (7.16) 
and it results in (primes dropped): 

% = -[^(^ + /)^ + ^vK mi. % = , = , 

= e-(^+^)^/2[(^/^)^a _ (2c)-i(6 + /)5° + 7(2c)-i(6 + /)5°] = nb'^. (7.17) 

In the new coordinates the metric is independent of t and y. After the coordinate trans- 
formation (7.16) and after solving the Killing equations for k'^-^ it assumes the form: 

^00 = \{b^ /)^e(^+/)^/.22, 501 = \{b^ De^'^'^^/'h,, 
902 = e^'+f>[-2c/{b + /) + 1(6 + f)Wh2 + ^(6 + fhhs], 
go3^l{b + f)e^'+f>/'h2s, 511 = /ill, 

= e('+^)"{-8c(6 + f)-^W + [2c/(6 + f)f + 2^Wh23 + l^h^^ + W^'/i22}, 

523 = e(''+^)^/2(W^^23 + 7/i33), 533 = /l33, (7.18) 

where the hij{z) are arbitrary functions, and W , not to be confused with the same symbol 
from sec. IV, is: 

W -.^-{b^ f )cx+ c. (7.19) 
2 



VIII. Case 2.1.2: c = 0. 

It follows from (7.1) that in this case: 



f = a^Q. (8.1) 
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Similarly as it happened in sec. VI, it turns out that this case is included in the case c 7^ 0. 
When the procedure of sec. VII is retraced with (8.1), the Killing fields that result at the 
stage corresponding to (7.14) are equivalent to (7.14), with 6 + / = 26, e playing the role 
of c and the roles of /c(2) and A;(3) interchanged. 

IX. Case 2.2: det A = 0. 

With A = det A = 0, the following follows from (2.8) and (2.6): 

bf — ce, f — —b, ai — a2 — 0, (9.1) 

and the matrix A is nilpotent, A^ — 0. Again the case c — has to be set aside for separate 
consideration, so we first follow: 

Case 2.2.1: c 7^ 0. 

The solutions of (2.5) - (2.6) are here: 

A2 = L2(i/)+ cxLsiy), = -{b/c)\2 + L^{y), (9.2) 

(f) = F{y) + cxP{y) + ^(a6 + cd)x^y, 

^ = -(6/c)0 + P{y) + c-\ab + cd)xy + c-^{-ay - 6C2 - cCg). (9.3) 

In considering eqs. (2.4c) and (2.3/), the case ab + cd = has to be considered separately, 
so now we follow: 

Case 2.2.1.1: ab + cd^ 0. 

Eq. (2.4c) is now a polynomial of second degree in x, and it implies the following 
equations: 

j = a, h = -d, {ab + cd){yF,y+F + C2)-c^PP,y={cg-a^)y. (9.4) 
The integral of the last equation is: 

F ^hab^cd)-\cPf/y+{cg-a^)y\+A{ab + cd)/y-C2. (9.5) 
Eq. (2.3Z) reduces to the single equation: 

-c2 (L3,^P + UP,y ) + (a6 + cd) {yL2,y + L2) = 0. (9.6) 

In solving (9.6) the case P = has to be considered separately. However, although some 
of the intermediate steps of the calculation depend on the assumption P 7^ 0, in the final 
formulae for the Kilhng fields P = is achieved by a transformation of the (1.13) set. The 
result is identical to the one obtained with P = from (9.6) on. 
With P 7^ 0, the solution of (9.6) is: 

L3 = {ab + cd)[yL2--i{ab + cd)]/{c^P). (9.7) 

We change the basis of the resulting Killing fields by taking k'^-^ = A;^-) + (a/c)A;"^) + (6/c)A;(^-), 
then carry out the transformation (1.13) with: 
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H ^cP/[y{ah + cd)], T,y= cL^/[y{ah + cd)], (9.8) 

and again change the basis to /c^g) = [cih + cd)"^[A;(^) — \{ah + cd)~^{cg — a^)/!;^*^)], A;""3) = 
c{ah + cd)~^k'^y The final KiUing fields are: 

^Td = '^T, ^-"2) = (2^/l/)'5o + (-^/l/' + ^72)5? - xy6^ + (7/l/)<53", 

^^3) = x5'^ - y5^. (9.9) 

The commutation relations are: 

k(2)] = [A';(2), A;(3)] = -%), [%), %)] = -%), (9.10) 

and they correspond to Bianchi type VIII. 

The Killing equations for imply that the metric is independent of x. Those for /c(3) 
are solved by: 

gii = Huy^, gi2 = Hu, gi3 = Hisy, 

5^22 = ^^22/y^ g23 = H23/y, ^33 = -f^33, (9.11) 

where the Hij{t, z) are arbitrary functions (the components g^a are as in (1.9), we are still 
in the Plebahski class in this case). In solving the Killing equations for k{2) it is useful to 
observe that A and 7 cannot vanish simultaneously: with ^ = 7 = 0, the Killing equations 
for k{2) imply that gi2 — g22 — 5'23 = 0, i.e. det{gap) — 0. With + 0, the following 
variables can be introduced: 

u = -ft- 2Az, v = 2At + -fz. (9. 12) 

The Killing equations for fc^^ have to be solved separately for Aj^O and for ^ = 0. When 
^ 7^ 0, we define 6 and A by: 

52 ^ 4^2 ^ ^2^ ^2 _ _2j[/s\ (9.13) 

and then the solution is: 

U :— hi2 sinh(2Av) + ki2 cosh(2Av), V :— /i23 sinh(Av) + k23 cosh(Av), 
Hn = (-2A)-'^^U + (7M)(-2^)-VV + hn, 
H12 = hi2 cosh(2A^;) + ku sinh(2AT;) + 7if23/(2^), 
Hi3 = (-2^)-VV + 7/133/(2^), H22 = (-2^)V2[/ + 2Ahn - ^^hs3/i2A), 

H23 = h23 cosh(At;) + A;23 sinh(Af ), H^s = h^s, (9.14) 

where the hij{u) are arbitrary functions. Eqs. (9.14) are adapted to the case A^ > 
(.4 < 0), but the solution for A^ < can be easily constructed from this one. 

When ^ = 0, necessarily 7 7^ and the Killing equations for /c^g) solved in the 
original variables {t, z) as follows: 

Hii = /i33^V(27)^ + /i23^V7^ + (/i22/7^ + /ii3/7)^^ + 2/ii2^/7 + hn, 
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Hi2 = h33Z^/{2-f) + 3/123-^7(27) + (^22/7 + hi3)z + hu, 

Hl3 = h33Z^/ (27) + h23z/j+ hi3, H22 = /isS-^^ + + /i22, 

-f^23 = h33Z + h23, H33 = /i33, (9.15) 

where the hij{t) are arbitrary functions. 

Eqs. (9.14) and (9.15) are very similar in form to (5.27) and (5.28), respectively. This 

suggests that Case 2.2.1.1 considered here may be included in Case 1.1.2.2 of sec. V as a 
limit (combined with a coordinate transformation). However, this author was not able to 
prove or disprove this hypothesis. 

We go back now to (9.3) and consider: 



X. Case 2.2.1.2: ab + cd = 0. 

Eqs. (9.2) still apply and eqs. (9.3) simphfy in the obvious way. The difference with 
sec. IX is that eqs. (2.4c) and (2.31) impose weaker conditions here. The case P — leads 
to the domain of Paper 2, so it may be assumed that P 7^ 0. Then, the consequences of 
(2.4c) are: 

h = bj/c, (10.1) 
cPP,y +(j - a)P +{g- aj/c)y = 0. (10.2) 
Eq. (10.2) defines P as a function of y. The consequence of eq. (2.3/) is: 

c(L3P),y+{j-a)L3^0, (10.3) 

which defines L3 once P is given. The functions F and L2 are still arbitrary at this point. 

In the resulting Killing fields we take k'^^-^ — k^^-^ + {a/c)k'(\^ + {b/c)k^2) instead of /c^^, 
and carry out the transformation (1.13) with: 

H = {F + C2)/{cP), T,y = {L2/C - HL3)/P. (10.4) 
The result of the change of basis and of (10.4) is equivalent to: 

F = -C2, L2 = a = 6 = d = 0, c=l, (10.5) 
and the Killing fields become: 

kti) = ^1. k^2) = 4P - yP,y K + xP,y - P^2 + ^L35t, 

k^^ = (P - yP,y + P,y + L36I (10.6) 

Note that with (10.5) fulfilled, the subcase g = implies P,y = const, and then /c^^ — 
P,y /c"^^ is spanned on and w", i.e. this case is in the domain of Paper 2. Hence, in 
what follows it will be assumed that g 0. For the same reason it will also be assumed 
that(P-yP,^)^0. 

The commutation relations (with (10.5) already taken into account) are: 

[fc(i), fc(2)] = %), [k(2), k(3)] = jk(3) + gkil), [%), = 0. (10.7) 
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The collection of Bianchi types contained in (10.7) is as follows. When: 

(a) g < the Bianchi type is VI/j with the free parameter = j/ — '^qY^'^'i and 
VIq when j = 0. 

(b) g > the Bianchi type is VII/j with the free parameter qb — j/ — j^)^/^; and 
VIIo when j — 0. 

(c) g = the Bianchi type is IV if j 7^ and II if j = 0. This last case belongs to 
the domain of Paper 2, as explained above, and will not be presented here. 

(d) g = with no condition on j, the Bianchi type is III. 

The Kilhng equations have to be solved separately for these three cases. 
The coordinates are adapted to and /c(3) by the following transformation: 

t' = t/{P - yP,y), x'^x- tP,y/{P - yP,y), y' = y, 

z' ^-Utl{P -yP,y)^z (10.8) 
In the new coordinates, with primes dropped: 

^fi) = ^"2) = - ^■^)^o - 9^^ - = ^0, 

u^^{P- yP,y)-\S^ - P,y - LsS^), = r^53^ 

^00 = - y'P,/ + + 2LsP,ygi3 + Ls^gss, 

901 = y(P - yP,y ) + P,y 911 + -^3«/l3, 
902 = P,y 9l2 + -f'35'23, 5^03 = P,y 9l3 + L^92.3. (10.9) 

The Killing equations for k{i) and A;(3) imply now that the metric is independent of t and 
X, while those for A;(2) are: 

-Pgii,y + 2y{P - yP,y ) + 2P,y gn + 2L^g^^ = 0, 

-P9l2,y + L-ig23 = 0, -P9l3,y + P:y 913 + -t'35'33 = 0, 

-Pg22,y - 2P,y g22 = 0, -Pg2s,y - P,y 523 = 0, -P^aa,, = 0. (10.10) 
The last three have the same solutions in all three cases: 

5-22 = /i22/P^, 5-23 = ^23/^, 5'33 = ^33, (10.11) 

where the hij{z) in the formulae (10.11) - (10.21) are arbitrary functions. The remaining 
equations in (10.10) have to be solved separately for each case. In each case it is useful to 
introduce the new variable Y{y) by: 

Y,y=l/P. (10.12) 

Then eq. (10.3) (with (10.5) taken into account) has the solution: 

La = (7/P)e-^'^. (10.13) 

In order to find P one has to differentiate (10.2) (again with (10.5) taken into account) by 
Y and use y,Y = P- The equation becomes: 
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P,YY+jP,Y+gP^O, (10.14) 

and y is then found from y — J PdY. The solutions of (10.10) and of (10.14) are now as 
follows: 

Ceise I: g < f/^. 

p = Me^i^ + ATe^^^, y = (M//xi)e^^^ + {N/i^2)e^'^, (10.15) 



where M and N are arbitrary constants and: 



1^1,2 = -j72 + £i,2(r/4 - 9r\ £i = 1, £2 = -1. (10.16) 

It can be assumed that MN ^ because in both the cases M = and N — we are 
back in the domain of Paper 2 (then P,y — const). The solutions of the remaining KiUing 
equations are here: 



911 = + hnP^ - 27/ii3e'^^^P/[M(/xi - fi^)] + -f^hs^e'^'"^ /[M{^i, - 1^2)?, 

912 = - 7/i23e'^^''/[M(A^i - ^^2)^], ^13 = hsP - ^hs^e'^''' /[M{pi, - pi^)]. (10.17) 

Case II: g > f/A. 
Here we define: 

D= (^^-JV4)^/^ f/ := Mcos(DF) + iVsin(DF), 

:=Msin(Dr) - Arcos(L>r), (10.18) 
where M and N are arbitrary constants, and then the solutions are: 

p = e-i^/^U, y = (4^2 + f)-'e-^^'^{ADV - 2jU), 
911 ^y' + hiiP' + 2^hse-^^UV/[D{M^ + N^)] + ^^hs^e-^'' /[D^M^ + N% 
912 = + ih2^V/[D{M^ + N^)Ul gis = h^sP + ^hsPV/[D{M^ + N^)U\. (10.19) 

Case III: g = j^A 0. 

With M and N being arbitrary constants, the solutions for P and y are here: 

P = {MY + N)e-^^l^, y = -2P/j - AMe'^^l'^/f. (10.20) 

We can assume M 7^ because with M = again P,j, = const. Then, the solutions of the 
Killing equations are as follows: 

g,r = + h,,P^ - 2^h,e-^^/^P/M + ^^h.-.e-^'' /M\ 
912 = hu - 7h23/[M{My + iV)], 513 = h^^P - ^h^se-^^/^M. (10.21) 

Case IV: gi = 0, no condition on j. 
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The conclusion that (10.5) can be achieved by a linear transformation of the basis of 
Killing fields is still valid. With a — g — and c — 1, eqs. (10.2) and (10.3) are solved by: 

P^-jy + M, L3 = 7, (10.22) 

and the Killing fields are then: 

(10.23) 

The metric that is a solution of the Killing equations for the above is: 

511 = {2M/j)y - {M/j f + {jy - Mfhn + 2{'y/j){jy - M)h^^ + (7/jf /i33, 

gi2 = hi2 + jj—zTm)' 9i3 = {jy - M)hi3 + (7/^)^33, 

g22^h22/{3y-M)\ g23^h23/{jy- M), gas^h^s, (10.24) 
where the hij are arbitrary functions of the argument 

T = -ft- Mz. (10.25) 
With this, case 2.2.1 is exhausted. We go back to (9.1) and consider: 

XI. Case 2.2.2 : c = 0. 

With A = det A = c = the only clement of the matrix A that may be nonzero is e. 
The solutions of (2.5) - (2.6) are here: 

A2 = L2(y), A3 = exL2 + L3(y), (11.1) 

(p^axy + F{y), 4^ = ^aex'^y + {dy + eC2 + eF)x + P{y). (11.2) 

Eq. (2.4c) now implies that either j = a or a = or e = 0. However, j = a leads to a 
result equivalent to (9.9) and a — leads to a result equivalent to (10.6). Hence: 

e = 0, (11.3) 
which makes the whole matrix ^4 = 0. The remaining implications of (2.4c) and (2.3Z) are: 

ah + dj - 0, y{dF,y -aP,y -g) = h{F + C2) + j{P + C3), 

y{dL2,y - aL^^y) = hL2 + JL3. (11.4) 
In solving these conditions, again several cases have to be considered separately. 

Case 2.2.2.1: a 7^ 0. 

Then: 

h — —dj/a, (11-5) 
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but in solving the second of (11-4) the case j — —a has to be set aside for separate 
consideration. We first follow: 

Case 2.2.2.1.1: j ^ -a. 

Then: 

P = -gy/{j + a)+ By-^''^ - C3 + (rf/a)(F + C2), 

L3 = -Ay-''^ + (rf/a)L2, (11.6) 
where A and B are arbitrary constants. In the resulting Killing fields we change the 

basis to A;(2) = a~"^^(2)' ^"(3) ^ ^fs) ~ (^/'^)^(2) + + "then carry out the 

transformation (1.13) with: 

H^{C2 + F)/{ay), T.^^Uj^ay). (11.7) 
The final Killing fields are: 

^fi) ~ ^(2) ~ ~ v^2i 

= y-^/lS(j7a + - Bj{ay)-'S^ - AS^]. (11.8) 
The commutation relations are: 

[k(i),k(2)] = [k(2), %)] = {j/a)k(3), [%), = 0, (11.9) 

and they correspond to Bianchi type Vlh with the free parameter as = {j — Ci)/{j + a). 
The coordinates will be adapted to and /c(3) after the transformation: 

t' = aty^/y[B{j + a)], x' = jt/[{j + a)y] + x, y' = y, z' = Aat + B{j + a)z. 

(11.10) 

After the transformation: 

kfi) = S^, k%) = -3t5^ + - ay5^, k^^ = ^o", 

u'^ = (j + a)-i[aS-iy^/"5- + {3/y)5'i] + ^a^g", = B{j + 0)7153", 

<7oo = {By-'/''f[l - {j/af + {j/afg^^/y^ + 2^j(j7« + l)^?i3/l/ + A^j + 0)^^33], 

^01 = %"^/"[(j7a + 1)1/ - {3la)giily - A{j + 0)^13], 

go2 = -By-^/''[{j/a)gi2/y + ^(j + 0)^23], ^03 = -By-^'''[{j / a)gi:,/y + A{j + a)(?33], 

(11.11) 

and the Killing equations imply: 

gu = 5'i2 = hi2, gi3 = ^i3y, 

5-22 = ^22/Z/^ 5'23 = ^23/y, 5-33 = ^33, (11-12) 

where the hij{z) are arbitrary functions. 
Case 2.2.2.1.2: j = -a. 
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The formulae for h and L3 simply follow from (11.5) and (11.6) while the second of 
1^11.4) now has a different integral: 



g/a)ylny+ By - + {d/a){F + C2). (11.13) 

(2) = ^ ^^(2)' ^(3) = ^(3) ~ {d/ci)k'^2) 



In the resulting Killing fields we change the basis to k'(2) = a ^kZ), k'/^^ = k%^ — {d/a)k?2) + 



{g/a — B)k"i^, and then carry out the transformation (1.13) with the same H and T as in 
(11.7). The final Killing fields k([-j and k^2) same as in (11.8), and: 

k^,^^y5^-\ny6'^-Ay5l (11.14) 
The commutation relations are: 

[^(1), ^(2)] = k(^i), [k(^2), ^(3)] = - fc(3), [fc(3), fc(i)] = 0, (11.15) 

and they correspond to Bianchi type IV. 

In order to adapt the coordinates to /c"^^ and /c^-j we now carry out the transformation: 

t'^t/y, x' ^{t/y)\ny + x, y'^y, z'^At + z, (11.16) 
and it leads to: 

^(1) ~ '^l"' ^(2) ~ ^'^O + (2^ ~ ~ y^2J ^(3) ~ '^O ' 



9oQ = - 2yMny + (lny)^5rii + 2Ay\nygi3 + {Ay f 933, 

901 = - In ygii - Ay 913, 902 = - In ygu - Ay 923, 

go3 = - \nygi3 - Aygss: (11-17) 

and the Killing equations lead to formulae for the other metric components that are iden- 
tical to (11.12). 

Case 2.2.2.2: a = 0. 

Then the first of (11.4) imphes that either d — 01 j — 0. However, if d = 0, then 
either the algebra of the Killing vectors necessarily becomes two-dimensional (this case is 
not considered here) or 

j = (11.18) 

follows anyway. Hence only (11.18) will be considered further. Then, with ^ d ^ h, 
a result equivalent to (11.8) foUows, and when d — h, a, result equivalent to Case 
2.2.2.1.2 follows. Therefore, new results are contained only in the case: 

d = Q. (11.19) 

Then, from (11.4): 

h{F + C2) + gy = 0. (11.20) 
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If /i 7^ 0, then (11.20) implies that the combination of the KiUing fields fc^-, + {g/h)k^i-^ is 
coUinear with w", and this case is in the domain of Paper 2. Hence, from (11.20): 

h = g = 0, (11.21) 

which leaves us with an Abelian algebra (Bianchi type I) and the functions F, P, L2 and 
L3 being all arbitrary. To avoid landing in the domain of Paper 2, we have to assume: 

F,y^O^Py, C2 + F-yF,y^0^C, + P-yP,y. (11.22) 

We will denote: 

WF:=C2 + F-yF,y, Wp -.^ Cs + P - yP,y , (11.23) 
and then the KiUing fields are: 

kfi) = 5?, kf2) = WpS^ + F,y 5? + A;f3) = WpS^ + P,y + L^S^ (11.24) 

It may also be assumed that: 

jC-.^Ls- L2Wp/Wf 7^ (11.25) 

because if £ = 0, then the Killing equations imply that either the symmetry group becomes 
two-dimensional or the metric is singular. With C ^ 0, the following transformation is 
permissible: 

t' = (Lgi - Wpz)/{jCWf), x' - {L2P,y -LsF,y )t/{LWF)+x- {P,y -F,y Wp/Wf)z/jC, 

y' = y, z' = -L2t/{CWF) + z/C, (11.26) 
after which the Killing fields become: 

kti) = S?, kf2) = Sl k^3)=Sl (11.27) 

The resulting metric depends only on y, but the Einstein equations are hopelessly compli- 
cated. 

This is the end of the collection of solutions of the KiUing equations. 
XII. Perspectives. 

The research for this series of papers was motivated by the desire to find a rotating 
(exact) perturbation of the Friedmann - Lcmaitre cosmological models. Several papers have 
been published whose authors found exact solutions of Einstein's equations with a rotating 
source (see a brief overview in sec. XIII), but all except one of them are either stationary 
from the beginning or become static in the limit of zero rotation. The one exception is 
the solution 2 by Stephani^ that has still nonzero expansion in the limit a; — > 0, but it is 
a rotating perturbation of a degenerate limit of the hyperbolically symmetric Kantowski - 



30 



Sachs solution^ (see Ref. 6 for more on this point) and cannot reproduce any Robertson - 
Walker geometry. 

There is no proof available that a rotating pertmbation of the F-L models can be 
spatially homogeneous. However, if it has any spatially homogeneous subcases and has a 
dust source, then the subcases must be contained among the metrics listed in this series 
of papers - just because this is a complete list of all hypersurface-homogeneous rotating 
dust metrics (and, over all three papers, all Bianchi types appeared on the list, some of 
them more than once). Hence, from now on, instead of testing various metric ansatzes for 
rotating dust by trial and error, one can choose an ansatz from a limited collection. It 
is well-known''' that the Robertson - Walker geometries have the following relation to the 
spatially homogeneous Bianchi-type geometries: 

The spatially flat {k — 0) R-W geometry is a common subset of Bianchi types I and 
VHo. 

The R-W geometry with negative spatial curvature (A; = —1) is a common subset of 
Bianchi types V and VH/j. 

The R-W geometry with positive spatial curvature {k — +1) is a subset of Bianchi type 
IX geometries. 

With this information, the following can be concluded: 

1. The metrics from Paper 1 do not contain any generalization of the F-L models; the 
Bianchi type I class contained there has timelike symmetry orbits and the velocity field is 
tangent to the orbits. 

2. The same is true for the Bianchi type 1 metrics from Paper 2 (cases 2.1.2.2 and 2.2) 
and from the present paper (eqs. (11.22) - (11.27)). This is in agreement with Theorem 3.1 
of King and EUis^ which says that no tilted Bianchi type 1 models exist (tilted means that 
the velocity field is not orthogonal to the symmetry orbits. Bianchi models with rotation 
obviously must be tilted). 

3. The Bianchi type V metrics from Paper 2 (cases 1.2.1 and 1.2.2.2) with suitably cho- 
sen parameters do have spacelike symmetry orbits and so may harbour some generalization 
of the k = —1 F-L model. 

4. The Bianchi type VH^ metrics of the present paper (eqs. (4.23) and (10.19)) may 
have spacelike symmetry orbits at least on some open subsets of the manifold, i.e. they 
may contain generalizations of the k = —1 F-L model. 

5. The same is true for the Bianchi type IX and type VHq metrics of the present paper 
(eqs. (5.23) with (5.27) or (5.28)); generalizations of the k = +1 and k = F-L models 
may be contained there. 

Hence, the cases listed in points 3, 4 and 5 are most promising from the point of view 
of cosmology. 

XIII. A brief overview of literature. 

Partly in order to justify the claim made in the first paragraph of sec. Xll, a brief 
overview of literature on solutions of Einstein's equations with rotating sources will be 
presented here. For the period up to 1973, the overview is based on a thorough survey of 
subject indexes to Physics Abstracts starting with the year 1915 (made in connection with 
Refs. 57 - 59). For the period 1973 - 1996 the survey was less thorough - I was looking for 
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only four keywords: "Bianchi", "homogeneous", "rotation" and "spatially homogeneous". 
In both searches the following sections of Physics Abstracts were surveyed: Cosmology, 
General Relativity, Gravitation and Space-Time Configurations (of course, references in 
the papers included in the survey were also checked). Ref. 9 is a more extended survey 
that covers the period up to 1973, and papers in which perturbative methods were used 
were also listed in it, these are omitted here. Vacuum solutions are omitted, too. 

As a rule, each paper is mentioned in one sentence, so this overview does not pretend 
to represent the contents of the papers; it is only meant to sort the papers by subjects and 
serve as an introductory guide through the literature. 

Lanczos^° found the historically earhest exact solution with rotating matter (although 
he may have been unaware of its rotation); it is a dust solution in which the velocity field 
and the rotation field are coUinear with Killing fields. Van Stockum^^ rediscovered the 
A = subcase of the Lanczos solution; Ref. 11 contains in addition important contributions 
to the techniques of solving the Einstein equations for stationary axisymmetric metrics 
(i.e. with a two-dimensional symmetry group) with a rigigly rotating source. Ref. 12 is 
another rediscovery of the Lanczos solution, and Refs. 13 and 14 contain discussions of its 
properties. 

Godel's solution^^ has a five-dimensional symmetry group whose orbit is the whole 
spacetime, and consequently the physical scalars in it are all constant. Its source is dust of 
zero expansion and shear, and constant matter density and rotation scalars. Rediscoveries 
of the Godel solution were published in Refs. 16 and 17 (see Ref. 18). Ref. 17 contains in 
addition a stationary axisymmetric perfect fiuid solution. 

Refs. 19 - 32 deal with properties of a metric form that is a modest generalization of 
the Godel solution and is known in the literature as the " Godel-type metric" . This notion 
was introduced in Ref. 19, and in Ref. 20 (which in fact preceded Ref. 19) it was shown 
that the only Godel-type metric with a perfect fluid source is the Godel solution itself. 
However, various Godel-type solutions with nonperfect fluid sources have been derived 
and investigated in Refs. 19 and 21 - 32. 

In Refs. 33 - 40 solutions of Einstein - Maxwell equations with a charged fluid source 
were discussed, some of them under the additional assumption that the Lorentz force is 
zero. The latter are fully within the scope of the formalism used in the present series 
of papers because the charged dust in them moves with zero acceleration. The relation 
of the results of these papers to those obtained here was described in Paper 1. Some of 
them are generalizations of the Lanczos^° and Godel^^ solutions. Those from Refs. 33 - 39 
are stationary cylindrically symmetric, i.e. are closely related to those from Paper 1, the 
one from Ref. 40 has a two-dimensional symmetry group. The solutions in Ref. 39 are 
coordinate transforms of those from Ref. 38. 

A generalization of the A = subcase of the Lanczos solution to a mixture of scalar 
field and dust was found by Santos and Mondaini^^. 

Other generahzations of the Godel solution were found in Refs. 42 - 50. Raval and 
Vaidya^^ found two solutions with anisotropic pressure, one of them nonstationary. Bray^^ 
found a collection of solutions with the rotating fiuid immersed in a magnetic field. Nov- 
ello and Rebougas^"^, Kitamura^^^^'' and Ray^^ found generalizations with heat fiow (the 
second of Kitamura's solutions presented in Ref. 45 is a coordinate transform of the Godel 
solution). Rebougas^^ found a generahzation with free electromagnetic field; see Ref. 51 
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and Paper 1 for an explanation of the relation between the solutions in Refs. 49 and 
37. Finally, Panov^'' found a generalization with two scalar fields, anisotropic fluid, null 
radiation and heat conduction. 

Other solutions with rotating charged matter source were obtained in Refs. 52 - 55. 
The Einstein - Maxwell equations for stationary axisymmetric charged dust were analysed 
by Bonnor^^. 

In addition, several solutions with a perfect fluid source (i.e. with nonconstant pressure) 
and with the same Bianchi type 1 symmetry as was considered in Paper 1 have been 
pubhshed. These include a family of solutions by this author^^~^° in which the velocity 
and the rotation fields are coUinear with the Killing fields. For the solutions considered in 
Refs. 57 and 59 the proportionality factor between -u;" and the Killing field is explicitly 
given. For the family of metrics from Ref. 60, the factor is an arbitrary function and the 
family is defined by a differential equation. The solutions by Davidson^^"^^ are explicit 
examples from the family of Ref. 60, possibly they are also coordinate transforms of 
members of the family from Ref. 57 defined by certain fixed values of parameters, but 
Refs. 61 and 62 do not contain sufficient information for precise identification. 

Nilsson and Uggla^^ did a qualitative analysis (using the theory of dynamical systems) 
of perfect fiuid solutions with the same symmetry that obey the linear barotropic equation 
of state. 

The list above includes papers that are related to Paper 1 of this series. Results related 
to those of Paper 2 are contained in Refs. 64 - 77. 

Among the solutions found by Ellis^^ there are some that directly belong to the col- 
lection of Paper 2, they are identified and described in Paper 2. Also within the scope of 
Paper 2 are the results of King^^ who investigated properties of the subcase /? = /113 = of 
Case 2.1.2.2 (of Paper 2) and provided a few examples of explicit solutions. Other explicit 
solutions in King's class were found by Maitra^^, Zimmerman^^ (this reference was not 
given in Paper 2) and Vishveshwara and Winicour^^. 

In Refs. 69 - 73 rotating dust solutions with four-dimensional symmetry groups were 
found; at least some of them have 3-dimcnsional subgroups and are within the domain 
of Paper 2. However, as mentioned in Paper 2, these papers do not contain sufficient 
information for a complete identification of all such subcases. The first three of the six 
solutions given in Ref. 70 are among the Case 2.1.2.2 metrics of Paper 2. 

Davidson'''^ found an example of solution with King's subclass of symmetry (stationary, 
cylindrically symmetric, differentially rotating) and with a perfect fiuid source obeying the 
linear barotropic equation of state. Two other stationary cylindrically symmetric perfect 
fiuid solutions were found by Garcia and Kramer^'', the first is differentially rotating and 
has the symmetry of King's subclass, the other is rigidly rotating and so has the Bianchi 
type I symmetry of the class considered in Paper 1. Nilsson and Uggla^^ analysed by the 
method of dynamical systems the Einstein equations for a metric with a Bianchi type II 
symmetry and a perfect fiuid source that obeys the linear barotropic equation of state. Its 
set of Killing fields is in the subcase A3 = of Case 1.1.2.2 of Paper 2. Stewart and Ellis''''' 
considered perfect fluid generalizations of the Ellis (dust) solutions from Ref. 64, they also 
considered sources with anisotropic pressure, viscosity and electric charge. 

The only paper directly relevant to the present Paper 3 is Ref. 3, see sec. V. The 
remaining part of the present section is a list of papers in which various problems connected 
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with rotating matter were discussed, but which are more remotely related to the present 
series of papers. 

Stephani"^ found a solution that is unique in one more respect in addition to that men- 
tioned in sec. XII: it is so far the only rotating matter (dust) solution with no symmetry. 
(Note: solution 3 in Ref. 4 is not a perfect fluid solution, contrary to the paper's state- 
ment). 

Stationary axisymmctric (i.e. with two-dimensional symmetry groups) perfect fluid 
solutions were found in Refs. 78 - 83. Wahlquist^* found a rigidly rotating solution with 
the equation of state e = —3p + const that, with specific values of two parameters, can 
describe the interior of a body with compact outer surface. Herlt^^ found a source of 
the NUT vacuum solution and Kramer^° found another solution, both are rigidly rotating 
(Kramer's solution was rediscovered by Patra and Roy^^, see Ref. 82). Another class of 
rigidly rotating stationary axisymmctric solutions was found by Herlt*^. 

Several stationary axisymmctric metrics were devised as nonperfect fluid sources of the 
Kerr solution. Results in this class that were published up to 1976 are reviewed in Ref. 
84. Later, a few more papers on this subject were published, but in those that are known 
to the present author the source is either a surface distribution of matter or an energy- 
momentum tensor that does not correspond to any identiflable kind of matter, hence they 
are not mentioned here. 

A spatially homogeneous solution of Bianchi type VIq with a rotating perfect fluid 
source was found by Rosquist®^; it has the equation of state e — 3p and nonzero expansion. 
Spatially homogeneous solutions with heat-conducting sources were found in Refs. 86 - 89. 
Two other hypersurface-homogeneous perfect fluid solutions were found by Wainwright^" , 
one has a 3-dimensional symmetry group of unidentifled Bianchi type, the other has a 
4-dimensional multiply transitive symmetry group. 

The remaining papers (Refs. 91 - 116) contain results obtained without flnding explicit 
solutions of Einstein's equations. This part of the survey is likely to be incomplete. 

Narlikar^^ proposed a metric form for a rotating dust model in which dust particles 
move on 3-cylinders (axial symmetry is not assumed) whose orthogonal sections are 2- 
dimensional surfaces of constant curvature. Winicour^^ reduced the Einstein equations 
with a stationary axisymmctric dust source to a sequence of integrations. This author-'-'^''- 
presented results of partial integration of the Einstein equations for a cylindrically sym- 
metric nonstationary perfect fluid. 

Bampi and Cianci^°^ investigated spacetimes with an Abelian 2-dimensional group of 
symmetries that has null orbits; the example of exact solution provided is a vacuum. Wils 
and van den Bergh^°^ showed that a stationary axisymmctric differentially rotating charged 
dust has cither a nonvanishing Lorentz force of a nonconstant ratio of charge density to 
mass density. 

In Refs. 105 - 110 general properties of rotating spatially homogeneous Bianchi type 
IX models were investigated without attempting to solve the Einstein equations. Of these, 
Ref. 109 gives a kinetic theory description of sources in such models and Ref. 110 gives a 
qualitative analysis of rotating mixmaster models. 

Refs. Ill - 113 contain general considerations about rotating matter models, and Refs. 
8 and 114 contain more general results, applicable also to nonrotating models, but having 
consequences for rotating models as well. CoUins^^^ considered properties of shearfree 
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rotating perfect fluids. Again CoUins^^^ reviewed arguments for the hypothesis that shear 
being zero imphes that either rotation or expansion is zero (apart from several specific 
examples, this is still unproven). Mason and Pooe^^^ investigated properties of the Lie 
derivative of the rotation vector along the velocity vector in rigidly rotating matter. King 
and Elhs^ investigated properties of tilted spatially homogeneous models and Nilsson and 
Uggla^^^ investigated hypersurface homogeneous and hypersurface self-similar perfect fluid 
models by the method of dynamical systems. 

Hawking^^^ and Collins and Hawking^^^ investigated limits set on the rotation param- 
eter in Bianchi-type models by the observations of the CMB radiation. 

Finally, there is a paper in the domain of the history of science. EUis^^^ described 
the influence of Godel's ideas presented in Refs. 15 and 105 on the development of several 
concepts and research programs in relativity, such as, among other things, the Bianchi-type 
models, singularity theorems and causal structure of spacetime. 

Acknowledgement Calculations for this paper were done with use of the algebraic 
computer program Ortocartan^^^"^^^. 

Note added in proof. The "Godel-type metric of Refs. 19 - 32 was discussed also in 
Ref. 120. Other rotating solutions of Einstein - Maxwell equations were discussed in Refs. 
121 - 122. I am grateful to M. J. Rebougas and A. Georgiou for this information. 

APPENDIX A 

The coordinate transformation (4.21) changes the metric tensor from the Plebahski 
form (1.9) to one in which the following relations hold (primes dropped, all the components 
displayed are expressed in the new coordinates): 

goo = e(^+^)^[sin^(L'x/2) - {V/Wf cos'{Dx/2)] + {VjWfg^^ 
+2(l^/W^)7e(''+-^)^/2[cos(L'x/2) + (l^/W^) sin(L'x/2)]^23 
+72e(^+^)^[cos(Da;/2) + iy jW) sm{Dx/2)fg^:i, 

901 = {V/W)gi2 + 'ye^'+^^^/^[cos{Dx/2) + {V/W) sm{Dx/2)]g^s, 

= e(''+-^)^cos(Dx/2)[sin(Dx/2) - (V/W) cos(Dx/2)] + {V/W)g22 
+7e(''+-^)^/2[cos(Dx/2) + {V/W) sin(Dx/2)]c/23, 

503 = {V/W)g2z + 7e('+^^"/'[cos(Da;/2) + {V/W) siYi{Dx/2)]g^^. 

APPENDIX B 
The equivalence of case 1.2 and case 1.1.1. 
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These are the main points of the reasoning from (2.8) on in the case c — 0. When c — 0, 
eqs. (2.7) - (2.8) imply: 

ai = b, a2 = f, (B.l) 
and b ^ f since A 7^ 0. The calculation requires checking several cases separately. 
Case 1.2.1: det A ^ 0. 

The formulae corresponding to (2.9) and (2.13) are then: 
A2 = Hy)e'^, A3 = [e/{b - f)]Hy)e'^ + L,{y)ef', <t> = F{y)e'^ - ay/b - C2, 

^ = [e/(6 - f)]F{y)e'^ + P{y)ef^ - [-ae/{bf) + d/f]y - C,. {B.2) 

In verifying (2.4c) and (2.3/) the cases F = and P = have to be considered 
separately. The results are as follows: 

When F — P and L2 — the group becomes 2-dimensional. 

When F = 7^ P and L2 7^ 0, the KiUing field [k^^^ + {a/b)kf^^] is coUinear with w'', 
and this case is in the domain of Paper 2. 

When P = 0, another linear combination of the Killing vectors with constant coefficients 
is coUinear with w". 

When Fy^Oy^P, formulae equivalent to (2.24) - (2.29) result. 

Case 1.2.2: det A = 0. 

Two different (but equivalent) subcases have to be considered here: 6 = and / = 
{by^f because of the assumption A 7^ 0). 

When f — 0, the functions A2, A3 and are given by the hmit / = of (B.2), while t/j 

is: 

^ = (e/6)F(y)e^^ + P{y) - aey/b'^ + {d - ae/b)xy - C3. {B.3) 

Eqs. (2.4c) and (2.3/) show then that there always exists a linear combination of the Killing 
fields with constant coefficients which is spanned on and w". 

When 6 = 0, the conclusion is the same, only the functions in (B.2) are different at the 
starting point. A2 and A3 are the hmits 6 = of those from (B.2), and: 

= F{y) + axy - C2, 
^ = -(e//)(0 + C2) + P{y)ef^ - {d/f + ae/f)y - C^. {BA) 
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CAPTION TO THE DIAGRAM 

The classes of metrics considered in the paper. Arrows point from more general classes 
to subclasses. The numbers at arrows are the case- numbers used in the text. The first 
entry in each rectangle is the property defining the case; all the symbols are introduced in 
cqs. (2.1) - (2.8). The subsequent entries give the following information: 1. The Bianchi 
type of the corresponding algebra (2.2); 2. The equation-numbers corresponding to the 
final result in the given case. No progress was made with the Einstein equations in any of 
the cases. Apart from case 1.1.2.2 for which a subcase was discussed in Ref. 3, none of the 
cases seem to have appeared in earlier literature (see sec. Xlll). 

The diagram does not show the links to the entries in the corresponding diagram in 
Paper 2; they are numerous and would obscure the drawing. 
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COMMUTATION RELATIONS 

- |A = 0| 



1.1.1.1,, 



- tti 
BVIh 
Eqs. (2.28) 
- (2.29) 



b + f^O 
BVIIh 
Eqs. (4.22) 
- (4.23) 











b + f=0 


BVIII and 




BIX, BVIII 


BVIq 




and BVIIq 


Eqs. (3.7) 
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ab + cd ifc 
BVIII 
Eqs. (9.9) - 
(9.15) 
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